SARS-CoV-2 infection elucidates unique features of pregnancy-specific immunity (preprint)

Pregnancy is a risk factor for increased severity of SARS-CoV-2 and other respiratory infections. The mechanisms underlying this risk have not been well-established, partly due to a limited understanding of how pregnancy shapes immune responses. To gain insight into the role of pregnancy in modulating immune responses at steady state and upon perturbation, we collected peripheral blood mononuclear cells (PBMC), plasma, and stool from 226 women, including 152 pregnant individuals (n = 96 with SARS-CoV-2 infection and n = 56 healthy controls) and 74 non-pregnant women (n = 55 with SARS-CoV-2 and n = 19 healthy controls). We found that SARS-CoV-2 infection was associated with altered T cell responses in pregnant compared to non-pregnant women. Differences included a lower percentage of memory T cells, a distinct clonal expansion of CD4-expressing CD8+ T cells, and the enhanced expression of T cell exhaustion markers, such as programmed cell death-1 (PD-1) and T cell immunoglobulin and mucin domain-3 (Tim-3), in pregnant women. We identified additional evidence of immune dysfunction in severely and critically ill pregnant women, including a lack of expected elevation in regulatory T cell (Treg) levels, diminished interferon responses, and profound suppression of monocyte function. Consistent with earlier data, we found maternal obesity was also associated with altered immune responses to SARS-CoV-2 infection, including enhanced production of inflammatory cytokines by T cells. Certain gut bacterial species were altered in pregnancy and upon SARS-CoV-2 infection in pregnant individuals compared to non-pregnant women. Shifts in cytokine and chemokine levels were also identified in the sera of pregnant individuals, most notably a robust increase of interleukin-27 (IL-27), a cytokine known to drive T cell exhaustion, in the pregnant uninfected control group compared to all non-pregnant groups. IL-27 levels were also significantly higher in uninfected pregnant controls compared to pregnant SARS-CoV-2-infected individuals. Using two different preclinical mouse models of inflammation-induced fetal demise and respiratory influenza viral infection, we found that enhanced IL-27 protects developing fetuses from maternal inflammation but renders adult female mice vulnerable to viral infection. These combined findings from human and murine studies reveal nuanced pregnancy-associated immune responses, suggesting mechanisms underlying the increased susceptibility of pregnant individuals to viral respiratory infections.

A single-nucleus and spatial transcriptomic atlas of the COVID-19 liver reveals topological, functional, and regenerative organ disruption in patients (preprint)

The molecular underpinnings of organ dysfunction in acute COVID-19 and its potential long-term sequelae are under intense investigation. To shed light on these in the context of liver function, we performed single-nucleus RNA-seq and spatial transcriptomic profiling of livers from 17 COVID-19 decedents. We identified hepatocytes positive for SARS-CoV-2 RNA with an expression phenotype resembling infected lung epithelial cells. Integrated analysis and comparisons with healthy controls revealed extensive changes in the cellular composition and expression states in COVID-19 liver, reflecting hepatocellular injury, ductular reaction, pathologic vascular expansion, and fibrogenesis. We also observed Kupffer cell proliferation and erythrocyte progenitors for the first time in a human liver single-cell atlas, resembling similar responses in liver injury in mice and in sepsis, respectively. Despite the absence of a clinical acute liver injury phenotype, endothelial cell composition was dramatically impacted in COVID-19, concomitantly with extensive alterations and profibrogenic activation of reactive cholangiocytes and mesenchymal cells. Our atlas provides novel insights into liver physiology and pathology in COVID-19 and forms a foundational resource for its investigation and understanding.

Longitudinal Characterization of Circulating Neutrophils Uncovers Distinct Phenotypes Associated with Disease Severity in Hospitalized COVID-19 Patients (preprint)

Multiple studies have identified an association between neutrophils and COVID-19 disease severity; however, the mechanistic basis of this association remains incompletely understood. Here we collected 781 longitudinal blood samples from 306 hospitalized COVID-19 + patients, 78 COVID-19 - acute respiratory distress syndrome patients, and 8 healthy controls, and performed bulk RNA-sequencing of enriched neutrophils, plasma proteomics, cfDNA measurements and high throughput antibody profiling assays to investigate the relationship between neutrophil states and disease severity or death. We identified dynamic switches between six distinct neutrophil subtypes using non-negative matrix factorization (NMF) clustering. At days 3 and 7 post-hospitalization, patients with severe disease had an enrichment of a granulocytic myeloid derived suppressor cell-like state gene expression signature, while non-severe patients with resolved disease were enriched for a progenitor-like immature neutrophil state signature. Severe disease was associated with gene sets related to neutrophil degranulation, neutrophil extracellular trap (NET) signatures, distinct metabolic signatures, and enhanced neutrophil activation and generation of reactive oxygen species (ROS). We found that the majority of patients had a transient interferon-stimulated gene signature upon presentation to the emergency department (ED) defined here as Day 0, regardless of disease severity, which persisted only in patients who subsequently died. Humoral responses were identified as potential drivers of neutrophil effector functions, as enhanced antibody-dependent neutrophil phagocytosis and reduced NETosis was associated with elevated SARS-CoV-2-specific IgG1-to-IgA1 ratios in plasma of severe patients who survived. In vitro experiments confirmed that while patient-derived IgG antibodies mostly drove neutrophil phagocytosis and ROS production in healthy donor neutrophils, patient-derived IgA antibodies induced a predominant NETosis response. Overall, our study demonstrates neutrophil dysregulation in severe COVID-19 and a potential role for IgA-dominant responses in driving neutrophil effector functions in severe disease and mortality.Funding:We acknowledge support from the Ragon Institute of MGH, MIT and Harvard, the Massachusetts Consortium on Pathogen Readiness (MassCPR), the NIH (3R37AI080289-11S1, R01AI146785, U19AI42790-01, U19AI135995-02, U19AI42790-01, 1U01CA260476 – 01, CIVIC75N93019C00052, T32 GM007592), the American Lung Association, and the MGH Executive Committee on Research. A.-C.V. acknowledges funding support from the COVID-19 Clinical Trials Pilot grant from the Executive Committee on Research at MGH, a COVID-19 Chan Zuckerberg Initiative grant (2020-216954), and funds from the Manton Foundation and the Klarman Family Foundation.Declaration of Interests: M.S.F receives funding from Bristol-Myers Squibb. G.A. is a founder of Seromyx Systems Inc. N.H. holds equity in Biontech and holds equity in and advises Danger Bio. Ethics Approval Statement: The study was approved by the Mass General Brigham Institutional Review Board under protocol 2017P001681, with an approval for a waiver of informed consent in compliance with the 45CFR 46, 2018 Common rule.

Longitudinal characterization of circulating neutrophils uncovers distinct phenotypes associated with disease severity in hospitalized COVID-19 patients (preprint)

Multiple studies have identified an association between neutrophils and COVID-19 disease severity; however, the mechanistic basis of this association remains incompletely understood. Here we collected 781 longitudinal blood samples from 306 hospitalized COVID-19+ patients, 78 COVID-19- acute respiratory distress syndrome patients, and 8 healthy controls, and performed bulk RNA-sequencing of enriched neutrophils, plasma proteomics, cfDNA measurements and high throughput antibody profiling assays to investigate the relationship between neutrophil states and disease severity or death. We identified dynamic switches between six distinct neutrophil subtypes using non-negative matrix factorization (NMF) clustering. At days 3 and 7 post-hospitalization, patients with severe disease had an enrichment of a granulocytic myeloid derived suppressor cell-like state gene expression signature, while non-severe patients with resolved disease were enriched for a progenitor-like immature neutrophil state signature. Severe disease was associated with gene sets related to neutrophil degranulation, neutrophil extracellular trap (NET) signatures, distinct metabolic signatures, and enhanced neutrophil activation and generation of reactive oxygen species (ROS). We found that the majority of patients had a transient interferon-stimulated gene signature upon presentation to the emergency department (ED) defined here as Day 0, regardless of disease severity, which persisted only in patients who subsequently died. Humoral responses were identified as potential drivers of neutrophil effector functions, as enhanced antibody-dependent neutrophil phagocytosis and reduced NETosis was associated with elevated SARS-CoV-2-specific IgG1-to-IgA1 ratios in plasma of severe patients who survived. In vitro experiments confirmed that while patient-derived IgG antibodies mostly drove neutrophil phagocytosis and ROS production in healthy donor neutrophils, patient-derived IgA antibodies induced a predominant NETosis response. Overall, our study demonstrates neutrophil dysregulation in severe COVID-19 and a potential role for IgA-dominant responses in driving neutrophil effector functions in severe disease and mortality.

SARS-CoV-2 Viremia is Associated with Distinct Proteomic Pathways and Predicts COVID-19 Outcomes (preprint)

BackgroundSevere Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) plasma viremia has been associated with severe disease and death in coronavirus disease 2019 (COVID-19) in small-scale cohort studies. The mechanisms behind this association remain elusive. MethodsWe evaluated the relationship between SARS-CoV-2 viremia, disease outcome, inflammatory and proteomic profiles in a cohort of COVID-19 emergency department participants. SARS-CoV-2 viral load was measured using qRT-PCR based platform. Proteomic data were generated with Proximity Extension Assay (PEA) using the Olink platform. ResultsThree hundred participants with nucleic acid test-confirmed COVID-19 were included in this study. Levels of plasma SARS-CoV-2 viremia at the time of presentation predicted adverse disease outcomes, with an adjusted odds ratio (aOR) of 10.6 (95% confidence interval [CI] 4.4, 25.5, P<0.001) for severe disease (mechanical ventilation and/or 28-day mortality) and aOR of 3.9 (95%CI 1.5, 10.1, P=0.006) for 28-day mortality. Proteomic analyses revealed prominent proteomic pathways associated with SARS-CoV-2 viremia, including upregulation of SARS-CoV-2 entry factors (ACE2, CTSL, FURIN), heightened markers of tissue damage to the lungs, gastrointestinal tract, endothelium/vasculature and alterations in coagulation pathways. ConclusionsThese results highlight the cascade of vascular and tissue damage associated with SARS-CoV-2 plasma viremia that underlies its ability to predict COVID-19 disease outcomes.

A single-cell and spatial atlas of autopsy tissues reveals pathology and cellular targets of SARS-CoV-2 (preprint)

The SARS-CoV-2 pandemic has caused over 1 million deaths globally, mostly due to acute lung injury and acute respiratory distress syndrome, or direct complications resulting in multiple-organ failures. Little is known about the host tissue immune and cellular responses associated with COVID-19 infection, symptoms, and lethality. To address this, we collected tissues from 11 organs during the clinical autopsy of 17 individuals who succumbed to COVID-19, resulting in a tissue bank of approximately 420 specimens. We generated comprehensive cellular maps capturing COVID-19 biology related to patients demise through single-cell and single-nucleus RNA-Seq of lung, kidney, liver and heart tissues, and further contextualized our findings through spatial RNA profiling of distinct lung regions. We developed a computational framework that incorporates removal of ambient RNA and automated cell type annotation to facilitate comparison with other healthy and diseased tissue atlases. In the lung, we uncovered significantly altered transcriptional programs within the epithelial, immune, and stromal compartments and cell intrinsic changes in multiple cell types relative to lung tissue from healthy controls. We observed evidence of: alveolar type 2 (AT2) differentiation replacing depleted alveolar type 1 (AT1) lung epithelial cells, as previously seen in fibrosis; a concomitant increase in myofibroblasts reflective of defective tissue repair; and, putative TP63+ intrapulmonary basal-like progenitor (IPBLP) cells, similar to cells identified in H1N1 influenza, that may serve as an emergency cellular reserve for severely damaged alveoli. Together, these findings suggest the activation and failure of multiple avenues for regeneration of the epithelium in these terminal lungs. SARS-CoV-2 RNA reads were enriched in lung mononuclear phagocytic cells and endothelial cells, and these cells expressed distinct host response transcriptional programs. We corroborated the compositional and transcriptional changes in lung tissue through spatial analysis of RNA profiles in situ and distinguished unique tissue host responses between regions with and without viral RNA, and in COVID-19 donor tissues relative to healthy lung. Finally, we analyzed genetic regions implicated in COVID-19 GWAS with transcriptomic data to implicate specific cell types and genes associated with disease severity. Overall, our COVID-19 cell atlas is a foundational dataset to better understand the biological impact of SARS-CoV-2 infection across the human body and empowers the identification of new therapeutic interventions and prevention strategies.

Codon arrangement modulates MHC-I peptides presentation: implications for a SARS-CoV-2 peptide-based vaccine (preprint)

Among various vaccination strategies, peptide-based vaccines appear as excellent candidates because they are cheap to produce, are highly stable and harbor low toxicity. However, predicting which MHC-I Associated Peptide (MAP) will ultimately reach cell surface remains challenging, due to high false discovery rates. Previously, we demonstrated that synonymous codon arrangement (usage and placement) is predictive of, and modulates MAP presentation. Here, we apply CAMAP (Codon Arrangement MAP Predictor), the artificial neural network we used to unveil the role of codon arrangement in MAP presentation, to predict SARS-CoV MAPs. We report that experimentally identified SARS-CoV-1 and SARS-CoV-2 MAPs are associated with significantly higher CAMAP scores. Based on CAMAP scores and binding affinity, we identified 48 non-overlapping MAP candidates for a peptide-based vaccine, ensuring coverage for a high proportion of HLA haplotypes in the US population (>78%) and SARS-CoV-2 strains (detected in >98% of SARS-CoV-2 strains present in the GISAID database). Finally, we built an interactive web portal (https://www.epitopes.world) where researchers can freely explore CAMAP predictions for SARS-CoV-1/2 viruses. Collectively, we present an analysis framework that can be generalizable to empower the rapid identification of virus-specific MAPs, including in the context of an emergent virus, to help accelerate target identification for peptide-based vaccine designs that could be critical in safely attaining group immunity in the context of a global pandemic.

Single cell profiling of COVID-19 patients: an international data resource from multiple tissues (preprint)

In late 2019 and through 2020, the COVID-19 pandemic swept the world, presenting both scientific and medical challenges associated with understanding and treating a previously unknown disease. To help address the need for great understanding of COVID-19, the scientific community mobilized and banded together rapidly to characterize SARS-CoV-2 infection, pathogenesis and its distinct disease trajectories. The urgency of COVID-19 provided a pressing use-case for leveraging relatively new tools, technologies, and nascent collaborative networks. Single-cell biology is one such example that has emerged over the last decade as a powerful approach that provides unprecedented resolution to the cellular and molecular underpinnings of biological processes. Early foundational work within the single-cell community, including the Human Cell Atlas, utilized published and unpublished data to characterize the putative target cells of SARS-CoV-2 sampled from diverse organs based on expression of the viral receptor ACE2 and associated entry factors TMPRSS2 and CTSL (Muus et al., 2020; Sungnak et al., 2020; Ziegler et al., 2020). This initial characterization of reference data provided an important foundation for framing infection and pathology in the airway as well as other organs. However, initial community analysis was limited to samples derived from uninfected donors and other previously-sampled disease indications. This report provides an overview of a single-cell data resource derived from samples from COVID-19 patients along with initial observations and guidance on data reuse and exploration.

Early immune response in mice immunized with a semi-split inactivated vaccine against SARS-CoV-2 containing S protein-free particles and subunit S protein (preprint)

The development of a vaccine against COVID-19 is a hot topic for many research laboratories all over the world. Our aim was to design a semi-split inactivated vaccine offering a wide range of multi-epitope determinants important for the immune system including not only the spike (S) protein but also the envelope, membrane and nucleocapsid proteins. We designed a semi-split vaccine prototype consisting of S protein-depleted viral particles and free S protein. Next, we investigated its immunogenic potential in BALB/c mice. The animals were immunized intradermally or intramuscularly with the dose adjusted with buffer or addition of aluminum hydroxide, respectively. The antibody response was evaluated by plasma analysis at 7 days after the first or second dose. The immune cell response was studied by flow cytometry analysis of splenocytes. The data showed a very early onset of both S protein-specific antibodies and virus-neutralizing antibodies at 90% inhibition regardless of the route of vaccine administration. However, significantly higher levels of neutralizing antibodies were detected in the intradermally (geometric mean titer - GMT of 7.8 {+/-} 1.4) than in the intramuscularly immunized mice (GMT of 6.2 {+/-} 1.5). In accordance with this, stimulation of cellular immunity by the semi-split vaccine was suggested by elevated levels of B and T lymphocyte subpopulations in the murine spleens. These responses were more predominant in the intradermally immunized mice compared with the intramuscular route of administration. The upward trend in the levels of plasmablasts, memory B cells, Th1 and Th2 lymphocytes, including follicular helper T cells, was confirmed even in mice receiving the vaccine intradermally at a dose of 0.5 g. We demonstrated that the semi-split vaccine is capable of eliciting both humoral and cellular immunity early after vaccination. Our prototype thus represents a promising step toward the development of an efficient anti-COVID-19 vaccine for human use.

RAPPID: a platform of ratiometric bioluminescent sensors for homogeneous immunoassays (preprint)

Heterogeneous immunoassays such as ELISA have become indispensable in modern bioanalysis, yet translation into easy-to-use point-of-care assays is hindered by their dependence on external calibration and multiple washing and incubation steps. Here, we introduce RAPPID (Ratiometric Plug-and-Play Immunodiagnostics), a "mix-and-measure" homogeneous immunoassay platform that combines highly specific antibody-based detection with a ratiometric bioluminescent readout that can be detected using a basic digital camera. The concept entails analyte-induced complementation of split NanoLuc luciferase fragments, photoconjugated to an antibody sandwich pair via protein G adapters. We also introduce the use of a calibrator luciferase that provides a robust ratiometric signal, allowing direct in-sample calibration and quantitative measurements in complex media such as blood plasma. We developed RAPPID sensors that allow low-picomolar detection of several protein biomarkers, anti-drug antibodies, therapeutic antibodies, and both SARS-CoV-2 spike protein and anti-SARS-CoV-2 antibodies. RAPPID combines ratiometric bioluminescent detection with antibody-based target recognition into an easy-to-implement standardized workflow, and therefore represents an attractive, fast, and low-cost alternative to traditional immunoassays, both in an academic setting and in clinical laboratories for point-of-care applications.

Plasma proteomics reveals tissue-specific cell death and mediators of cell-cell interactions in severe COVID-19 patients (preprint)

COVID-19 has caused over 1 million deaths globally, yet the cellular mechanisms underlying severe disease remain poorly understood. By analyzing several thousand plasma proteins in 306 COVID-19 patients and 78 symptomatic controls over serial timepoints using two complementary approaches, we uncover COVID-19 host immune and non-immune proteins not previously linked to this disease. Integration of plasma proteomics with nine published scRNAseq datasets shows that SARS-CoV-2 infection upregulates monocyte/macrophage, plasmablast, and T cell effector proteins. By comparing patients who died to severely ill patients who survived, we identify dynamic immunomodulatory and tissue-associated proteins associated with survival, providing insights into which host responses are beneficial and which are detrimental to survival. We identify intracellular death signatures from specific tissues and cell types, and by associating these with angiotensin converting enzyme 2 (ACE2) expression, we map tissue damage associated with severe disease and propose which damage results from direct viral infection rather than from indirect effects of illness. We find that disease severity in lung tissue is driven by myeloid cell phenotypes and cell-cell interactions with lung epithelial cells and T cells. Based on these results, we propose a model of immune and epithelial cell interactions that drive cell-type specific and tissue-specific damage in severe COVID-19.

Looking for pathways related to COVID-19 phenotypes: Confirmation of pathogenic mechanisms by SARS-CoV-2 - Host interactome. (preprint)

In the last months, many studies have clearly described several mechanisms of SARS-CoV-2 infection at cell and tissue level. Host conditions and comorbidities were identified as risk factors for severe and fatal disease courses, but the mechanisms of interaction between host and SARS-CoV-2 determining the grade of COVID-19 severity, are still unknown. We provide a network analysis on protein-protein interactions (PPI) between viral and host proteins to better identify host biological responses, induced by both whole proteome of SARS-CoV-2 and specific viral proteins. A host-virus interactome was inferred on published PPI, using an explorative algorithm (Random Walk with Restart) triggered by all the 28 proteins of SARS-CoV-2, or each single viral protein one-by-one. The functional analysis for all proteins, linked to many aspects of COVID-19 pathogenesis, allows to identify the subcellular districts, where SARS-CoV-2 proteins seem to be distributed, while in each interactome built around one single viral protein, a different response was described, underlining as ORF8 and ORF3a modulated cardiovascular diseases and pro-inflammatory pathways, respectively. Finally, an explorative network-based approach was applied to Bradykinin Storm, highlighting a possible direct action of ORF3a and NS7b to enhancing this condition. This network-based model for SARS-CoV-2 infection could be a framework for pathogenic evaluation of specific clinical outcomes. We identified possible host responses induced by specific proteins of SARS-CoV-2, underlining the important role of specific viral accessory proteins in pathogenic phenotypes of severe COVID-19 patients.

Induction of a regulatory myeloid program in bacterial sepsis and severe COVID-19 (preprint)

A recent estimate suggests that one in five deaths globally are associated with sepsis. To date, no targeted treatment is available for this syndrome, likely due to substantial patient heterogeneity and our lack of insight into sepsis immunopathology. These issues are highlighted by the current COVID-19 pandemic, wherein many clinical manifestations of severe SARS-CoV-2 infection parallel bacterial sepsis. We previously reported an expanded CD14+ monocyte state, MS1, in patients with bacterial sepsis or non-infectious critical illness, and validated its expansion in sepsis across thousands of patients using public transcriptomic data. Despite its marked expansion in the circulation of bacterial sepsis patients, its relevance to viral sepsis and association with disease outcomes have not been examined. In addition, the ontogeny and function of this monocyte state remain poorly characterized. Using public transcriptomic data, we show that the expression of the MS1 program is associated with sepsis mortality and is up-regulated in monocytes from patients with severe COVID-19. We found that blood plasma from bacterial sepsis or COVID-19 patients with severe disease induces emergency myelopoiesis and expression of the MS1 program, which are dependent on the cytokines IL-6 and IL-10. Finally, we demonstrate that MS1 cells are broadly immunosuppressive, similar to monocytic myeloid-derived suppressor cells (MDSCs), and have decreased responsiveness to stimulation. Our findings highlight the utility of regulatory myeloid cells in sepsis prognosis, and the role of systemic cytokines in inducing emergency myelopoiesis during severe bacterial and SARS-CoV-2 infections.